Coherent dynamics of domain formation in the Bose Ferromagnet 
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We present a theory to describe domain formation observed very recently in a quenched 87 Rb 
gas, a typical ferromagnetic spinor Bose system. An overlap factor is introduced to characterize 
the symmetry breaking of Mf — ±1 components for the F = 1 ferromagnetic condensate. We 
demonstrate that the domain formation is a co-effect of the quantum coherence and the thermal 
relaxation. A thermally enhanced quantum-oscillation is observed during the dynamical process 
of the domain formation. And the spatial separation of domains leads to significant decay of the 
Mf — component fraction in an initial Mf = condensate. 

PACS numbers: 05.30.Jp, 03.75.Kk, 03.75.Mn 



Very recently, the Berkeley group observed sponta- 
neous symmetry breaking in 87 Rb spinor condensates [l[ . 
Ferromagnetic domains and domain walls were clearly 
shown using an in-situ phase-contrast imaging. This ap- 
pears the first image of the domain structure in a Bose 
ferromagnet. Although ferromagnetism has been inten- 
sively studied in the context of condensed matter physics 
and is regarded as one of the best understood phenom- 
ena in nature Q, the description of ferromagnetism is 
not yet complete. The conventional ferromagnets being 
considered are usually comprised of either classical parti- 
cles (insulating ferromagnets) or fermions (itinerant fer- 
romagnets) while Bose systems are seldom touched 
The realization of cold spinor 87 Rb gases a typical 
ferromagnetic Bose system, has provided an opportunity 
to study Bose ferromagnets and thus opens up a way to 
a comprehensive understanding of ferromagnetism in all 
kinds of condensed matters. 

The ferromagnetic spinor Bose gas has attracted nu- 
merous theoretical interests 0, 0, HI, 0, EH U, 12, 13, 



FJ, |15|. On one side, researchers expect that it will show 



some general properties as conventional ferromagnets do. 
Ho Ohmi and Machida Q pointed out that this sys- 
tem has a spontaneous symmetry-broken ground state 
and a normal spin-wave excitations spectrum at small 
wave vector k, u> 3 = c s k 2 . On the other side, researchers 
aim at exploring distinct features of the system. Studies 
on thermodynamics and phase transitions have revealed 
that the ferromagnetic spinor Bose gas displays a quite 
surprising phase diagram. Its Curie point can be larger 
by magnitudes than the energy scale of the ferromag- 
netic interaction between bosons, and never below the 
Bose-Einstein condensation point 0, HI]. It means that 
once the Bose gas condenses, it is already spontaneously 
magnetized. 

A conventional ferromagnet usually has some domain 
structure below the Curie point, as illustrated in Fig. la. 
But whether it is true for a Bose ferromagnet is still some- 



what controversial. One even questions whether there 
exists a Curie point in ferromagntic Bose gases 0, [l(| , 
as the cold atomic gas under experimental conditions is 
usually not in the thermodynamic limit while the phase 
diagram mentioned above is derived from the thermal- 
equilibrium grand canonical ensemble 0, Hj]. Never- 
theless, a number of theoretical works have discussed 
the possibility of the domain formation 0, 14 , EH • 
Within the mean-field theory, Zhang et al. found out 
that the ferromagnetic condensate has a dynamical in- 
stability leading to spontaneous domain formation in an 
initially magnetized Bose gas [l3j]. Moreover, Mur-Petit 
et al. showed that a multi-spin-domain structure mani- 
fests in 87 Rb condensates at finite temperatures The 
Berkeley experiment confirms that the Bose ferromagnet 
can indeed form domain structures at least under certain 
conditions [l[. Then new questions come, is the process 
of the domain formation similar to that inside a conven- 
tional ferromagnet, and how does the domain formation 
affect spin dynamics? We attempt to answer these ques- 
tions in the present letter. 

We start with the Hamiltonian for anF=l system in 
the following form Q 



H = d 



+ J^lMa + f ' Fa'b^b^b] ,(1) 



where i/j a (r) is the field annihilation operator for an atom 
in state Mp = a at point r, /i is the chemical potential 
and U is the trapping potential, go and gi are the spin- 
independent and spin-dependent mean-field interaction, 
respectively. For macroscopically occupied Bose systems 
it is common to replace the field annihilation operator 
for the ath spin component by its expectation value, i.e. 
ip a (r, t) = {ip a {Y, t)), which for spinor condensates is con- 
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FIG. 1: Schematic domain structure inside a Ferromagnet. 
(a) The bulk material shows no magnetism because domains 
are randomly oriented, (b) The multi-domain configuration is 
often simplified to a two-domain structure with opposite sign 
of magnetization for theoretical convenience. The gray region 
denotes the domain wall, (c) Two-domain structure for a fer- 
romagnetic spinor Bose-Einstein condensate. Dotted, dashed 
and solid lines represent the normalized particle distributions 
of Mf = 1, Mf = — 1 and Mf = bosons respectively, (d) In 
a homogeneous system, the normalized particle distribution 
is taken as a constant in each domain; for Mf = bosons it 
remains a constant through the whole system. 



veniently expressed as 



<p a (r,t) = vW^e*^'' 1 



(2) 



Here N a (t) is the number of condensed particles, 
7? a (r, t) denotes the normalized particle distribution with 
j d 3 rrj a = 1 and <f> a {r,t) the phase. The multiplier 
of N a and r\ a refers to the condensed particle density 
of the a-component, and the total particle density is 
p(r,t) — N a (t)rj a (r,t). We suppose that the particle 
distribution can be different for distinct components in 
our approach. Neglecting excitations, the Hamiltonian is 
simplified to H = Hq + H u + H s , with 



Ha 



H„ = 



d 3 r 



d 3 r 
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(J-)N a rj a + —N a r] a N a ,ri a , 



d 3 rf [Nlrj* 



N 2 n 2 _ 

+2N r]o(N + r] + + AT_r ? _) 



2N + N-7] + r]- 



4N ^N+N-rio ^qIcos6 



(3) 



where 9 — <f>+ + <p- — 20o is the relative phase. 

In case that the ground state of the condensate is 
symmetry-broken, certain magnetic domain structure is 
formed spontaneously. Within each domain the atomic 
magnetic moments are aligned in a preferential direc- 
tion. A direct consequence of domain formation is that 
the Mp = 1 and Mp = —1 components are spa- 
tially separated, as portrayed in Fig. lc. The integral 



J d 3 r?7 + ?7_ measures the extent of overlap between the 
two, which is called the overlap factor hereinafter. In 
general, several "overlap factors" should be introduced 
to Eqs.(3), e.g., ao± = V J d 3 r-q±, a\ = V J d 3 rr] + r]-, 
ct2± = V J d 3 rr]Qrj± and = V J d 3 rr]oy/ri + r/_, where V 
is the volume of the system. Treating the relative phase 
as a spatially inde pen dent constant as previous theory 
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the term of H s in Eqs. (3) is 



rewritten as 



H„ = 



[a 0+ Nl + a -N 2 _ - 2 ai N+N_ 
+2N (a 2+ N + +a 2 -N-) 



(4) 



According to their definition, these overlap factors are 
not totally independent from each other and the number 
of independent ones can be further reduced. For simplic- 
ity, we consider a homogeneous spinor Bose gas with a 
two-domain structure, as shown in Fig. Id, and neglect 
the domain wall. In this case, we derive, after some in- 
tegration and algebraic manipulation, that there is only 
one independent overlap factor and the above equation 
is reduced to 



H s = ^ [(2 - a)(Nl +Nl)- 2aN+N^ 



+ 2N {N+ + N_) + A^Nny/N+N-cose , (5) 

where the reduced overlap factor is just given by a — 
V J d 3 TT] + ri_ . Our model allows the ev-factor to vary from 
one to zero, corresponding to the case that the two com- 
ponents are from thoroughly mixed to completely sep- 
arated. For a homogeneous system, the gradient term 
Hq can be doped. The spin-irrelevant term H u remains 
a constant since the total density distribution is ha rdly 
responsive to the evolution of domain structure 0, (Hf, 
p(r, t) w p(r, 0). Thus the dynamics of domain formation 
is only determined by the Hamiltonian expressed in Eq. 

The Berkeley experiment considered a pure spinor con- 
densate initially prepared in the unmagnetized state. It 
is important to emphasize that the total spin is con- 
served in an atomic quantum gas under experimental 
conditions ill 17 , IH ■ Therefore the particle numbers of 
Mp = 1 and —1 component are always equal, N+ = N-. 
Such a system is described by the Hamiltonian 



H s = -(1 - a)(l - n ) 2 - 2n (l - n ) [l + Vacos6»] .(6) 

Here H s = H„/(N\%\) with total particle number N = 
N + + N- + N \ n a = N a /N is the fraction of Mp = a 
component, no and form a pair of conjugate variables 
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FIG. 2: Contour plot of the energy surface in the 0-no plane 
with the overlap factor a = 0.5 (a) and a — 0.1 (b). (c) 
shows the fraction of Mf = atoms at the minima, no, as a 
function of a. 




FIG. 3: Dynamics of domain formation for Bose systems with 
the overlap factor a = 0.1 and Tr = 20, 10, and 2 from top 
to bottom. The left column shows the particle fractions of 
Mf = (no, solid lines) and Mf = ±1 (n±, dotted lines) 
components. The right column shows the magnetization of 
magnetic domains, m. The initial state configuration is (n+ = 
0.005, n = 0.99, n_ = 0.005). 



-n 



and their equations of motion are given by 

d_ 

at 

-2(1 - 2n ) [1 + \/acos6 



— 2 v / cmo(l — no)sin#, 
2(l-a)(l-no) 



(7a) 



(7b) 



The population dynamics depicts the process of the do- 
main formation; and the magnetization of magnetic do- 
mains is defined as m = N + r/ + — N—i]-. 

A number of theoretical works have investigated the 
spin dynamics of the 87 Rb spinor condensate II], II, HI • 
We notice that previous theories usually treated the three 
components being mixed as they share the same spatial 
wave function, which is known as the single-mode approx- 
imation. Therefore the domain structure is smeared out 
and thus the ferromagnetic feature was not sufficiently 
dealt with. Those theories correspond the special case of 
a = 1 in this letter. 

The dynamical behaviors of Eqs. (J7J) can be visual- 
ized by the phase-space portrait with constant energy 
lines. Figure 2a and 2b plot the contour lines of en- 
ergy with the overlap factor a — 0.5 and 0.1, respec- 
tively In case of a — 1. there are energy minima along 
no = 1/2 at 9 = 2nir |ll|, El El and all the contour 
lines are closed loops around those minima. As domains 



build up, the value of a drops down and two apparent 
changes take place, seen in Fig. 2a and 2b. (i) There 
exist two distinct regimes in the phase space diagram. 
The newly appeared regime lies in the upper region of 
the figures, which consists of a set of open curves. Each 
line corresponds to a rotation type of solution, in which 
the relative phase 9 is "running" with the time. The 
closed orbits, lying in the lower region, represent libra- 
tion type of solutions, in which 6 oscillates around the 
minimum, (ii) The positions of those minima move to- 
wards to smaller values of no. The minima points are 
connected to the ground state of the system. Figure 
2c shows the Mp = particle fraction at the minima, 
fiQ. The less a is, the smaller no is. It means that 
the domain formation tends to reduce the number of 
Mf = atoms. This point has also been affirmed by 
Mur-Petit et al. who obtained a state with equipartition 
in populations, (n+ ~ 1/3, no rs 1/3, n_ w 1/3), from 
a starting state (0.005,0.99,0.005) pj. We derive that 
no = n + = n_ = 1/3 when a drops to 0.25. In the limit 
case of a = 0, no = and n + = n_ = 0.5. 

The state of the quenched Mp = condensate is 
viewed as a point lying in the upper region of the phase 
space diagram (Fig. 2a or 2c), when Eqs. ([7]) yield a self- 
trapping solution [1 9(| . This motion reflects the quantum 
mechanical nature of the Bose-Einstein condensate. The 
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oscillating amplitude of no is very small and the resulting 
magnetization m is so low in magnitude that it is hard to 
be probed experimentally during this stage 0. The self- 
trapping effect prevents the growth of magnetic domains. 
This case is similar to the classical Larmor precession of 
a spin around magnetic field: it is rotating all the time, 
but the spin direction can be never parallel to the field 
without energy dissipation. 

Then one has to take into consideration the effect of 
thermal agitation, which can change the energy of the 
system and drive the system into thermal equilibrium. 
Therefore, if no departs from its thermal equilibrium 
value no > h will relaxes to uq exponentially with a charac- 
teristic time scale Tr, called the relaxation time. Assum- 
ing that the relaxation velocity is proportional to no— no, 
dfio/dt oc nn_— no, Eq. (fTaj) can be replaced with the fol- 
lowing one [20| | . 

^-n = -2Vano(l - n o )sin0 + "° "° . (8) 
at i R 

T R scales qualitatively the thermal dissipation rate. 
Longer T R denotes weaker thermal agitation. 

Figure 3 displays the population of Mp = and ±1 
components, as well as the magnetization m. The evo- 
lution of no and m reflects the dynamical process of the 
domain formation. As shown in Fig. 3a and 3b, no de- 
creases os dilatorily with time driven by the the thermal 
agitation, and meantime m arises. A very interesting re- 
sult is that the oscillation amplitude increases as the sys- 
tem relaxes. Generally, the thermal agitation suppresses 
the macroscopic quantum coherence, and thus tends to 
kill oscillations of the population, while here we see that 
the oscillation is enhanced. Furthermore, the magnetiza- 
tion persists in oscillating for a very long period of time 
after the amplitude reaches its maximum. This result is 
qualitatively consistent with the Berkeley group's obser- 
vation of the unstable magnetization mode[l[. 

If the thermal agitation gets stronger, the oscillation 
will be enhanced first, and then suppressed, as shown in 
Fig. 3c and 3d. Correspondingly, we divide the whole 
process into two periods with respect to the domain for- 
mation, the growing period and the stabilizing period. In 
the latter period, n and m oscillate around their thermal 
equilibrium values and the amplitudes decrease gradu- 
ally, then we have a stable domain structure eventually. 
If sketching the solution in the phase space diagram, one 
can find that the growing period is represented by the 
trajectory in the libration regime, and solutions for the 
stabilizing period lie in the rotation regime. Based on 
this understanding, Fig. 3a and 3b show only the grow- 
ing period. Given the thermal agitation strong enough, 
the quantum mechanical feature will be smeared out in 
both periods, as Fig. 3e and 3f indicate. 

According to the above discussions, the present model 
can qualitatively describe the dynamic process of the do- 
main formation in a quenched Mp — condensate. Sig- 



nificantly, we show that the spatial separation of mag- 
netic domains brings about much nontrivial effects on the 
spin dynamics of the ferromagnetic condensate. To get a 
quantitative description, more local details of the particle 
distribution and the relative phase should be considered. 

In conclusion, we have investigated the dynamics of do- 
main formation in a ferromagnetic spinor Bose-Einstein 
condensate, taking into account of the symmetry- 
breaking of the Mp — 1 and —1 components. Magnetic 
domains develop with the separation of Mp = ±1 compo- 
nents. Our results suggest that the Mp = component 
in the condensate can significantly decay to a very small 
value, far less than 1/2 as previous theories predicted. 
The domain structure is formed and stabilized with the 
help of the thermal dissipation. A thermally enhanced 
quantum-oscillation is observed during the process. 
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